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'i> '. Abstract 

\o . 

■ The HERA data on the proton structure function, F2{x, Q^), at very small x and show 
^ ■ the dramatic departure of the logarithmic slope, (9^2 /(9 log Q^, from theoretical predictions 

! based on the DGLAP evolution. We show that the running BFKL approach provides the 
^ I quantitative explanation for the observed x and/or -dependence of dF2/dlogQ^. 
^ ■ 

O . Caldwell's presentation of the HERA data in terms of the logarithmic derivative dF2 / d log 

for the proton structure function (SF) F2{x, Q^) exhibits the turn-over of the slope towards 
small X and/or up to currently attainable x ~ 10~^ and ~ 0.1 GeV^ @]. The 
^ ■ DGLAP-evolution with GRV input predicts a steady increase of the derivative 

H : 

. . . aj^DGLAP 

^^^-^^asmG^^'^^^{x,Q') (1) 
o log Q'^ 

with 1/x, due to the growth of the gluon structure function G^'^^^^{x, Q^) = xg^^^^^ {x, Q"^). 
A slight systematic discrepancy of the DGLAP analysis with small-x data on dF2/ dlogQ"^ 
was found also in 

The turn-over point located at x ~ 5 ■ 10^^ and ~ 5GeV^, in a commonly believed 
legitimate pQCD domain. So, the phenomenon occurs on the interface between "soft" and 
"hard" physics. Its explanation within the color dipole approach is based on two observa- 
tionsQ: 

i) specific smallness of the log Q^-derivative of sub-leading terms of the BFKL-Regge expan- 
sion for F2 at the turn-over point, which is due to the nodal structure of the running BFKL 
eigen-SF's; 

ii) significant contribution to the small-Q^ proton SF coming from the non-perturbative 
^The preliminary results have been reported at the DIS'98 Workshop 
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component of the BFKL pomeron. 



The s-channel approach to the BFKL equation was developed in terms of the color 
dipole cross section a{x,r) |^ (hereafter r is the color dipole moment). The positive 
feature of the color dipole picture, to be referred to as the running BFKL approach, is 
consistent incorporating the two crucial properties of QCD: i) asymptotic freedom (AF), 
i.e., the running QCD coupling as{r) and, ii) the finite propagation radius Rc of perturbative 
gluons. 

The BFKL equation for the interaction cross section a{x, r) of the color dipole r with 
the target reads 



d\og{l/x) 



3 
8^ 



d pi \S{pi) - S{p2)\ [a{x,pi) + a{x, P2) - cr{x,r)] . (2) 



Here the kernel /C is related to the wave function squared of the color-singlet qqg state with 
the Weizsacker- Williams (WW) soft gluon. The quantity 

^(p) = -9s{p)^pKo{pGp) = 9s{p)PgKi{pgP)p/p, (3) 

where Rc = l/pc and K^{x) is the modified Bessel function, describes a Yukawa screened 
transverse chromoelectric field of the relativistic quark and \S{pi) — S{p2)\'^ describes the 
flux (the modulus of the Poynting vector) of WW gluons in the qqg state in which r is 
the q-q separation and pi^2 Sive the q-g and q-g separations in the two-dimensional impact 
parameter plane. Our numerical results are for the Yukawa screening radius Rc = 0.27 fm. 
The recent fits to the lattice QCD data on the field strength correlators suggest similar Rc 
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The asymptotic freedom of QCD uniquely prescribes the chromoelectric field be com- 



puted with the running QCD charge gs{r) = yAnasir) taken at the shortest relevant 
distance, Ri = min{r, pj} in the qqg system. Although, the so introduced running coupling 
does not necessarily exhaust all NLO effects, it correctly describes the crucial enhancement 
of long distance, and suppression of short distance, effects by AF. 

The properties of the running color dipole BFKL equation responsible for the observed 
dependence of dF2/dlogQ^ are as follows jri], |12|. The spectrum of the running BFKL 
equation is a series of moving poles in the complex j-plane with eigen-functions 

anix,r) = cr„(r)exp [A„ log(l/a;)] (4) 

being a solution of 

/C (g) (T„ = A„cr„(r). (5) 

The leading eigen-function ao{r) is node free. The sub-leading o"„(r) has n nodes. The 
intercepts A„ closely, to better than 10%, follow the law A„ = Ao/(n -|- 1) suggested earlier 
by Lipatov |jl3|. The intercept of the leading pole trajectory, with the above specific choice 
of -Rc, is Aq = Aip = 0.4. The sub-leading eigen-functions o"„ ||TT], O] are very close to 



Lipatov's quasi-classical solutions [0 for n ^ 1. For our specific choice of the infrared 
regulator, Rc, the node of cri(r) is located at r = ri ~ 0.05 — 0.06 fm, for larger n the first 
node moves to a somewhat larger r ~ 0.1 fm. 
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The color dipole factorization [T^ in conjunction with the exphcit form of the qq hght- 

, relates the dipole cross sections o"„(r) with the eigen-SF, 



cone wave function, '^qq[z,r 



fn{Q' 



q=u,d,c,s 



The BFKL-Regge expansion 

a{x,r) = cro{r){xo/x)'^° + ai{r){xQ/x)^^ + o-2{r){xo/x)'^'' + .. 
gives the BFKL-Regge expansion for the SF 

F2{x,Q') = Y.UQ'){xo/^. 



'X 



(6) 



(7) 



(8) 



The remarkable finding of ^S], |Tl], |T2| is a good description of the HERA data on the 
proton SF starting with the Born two-gluon cross section (TBir) as a boundary condition 
for the running BFKL equation at Xq = 0.03. With such a boundary condition, which 
could well be excessively restrictive, the expansion (0) fixes uniquely the normalization of 
the eigen-FS's. 

The Bjorken variable x = Q"^ /2mpV is commonly being used for the presentation of the 
experimental data even at < m^, way beyond the kinematical region ^ it has 
originally been devised for. At small Q^, the relevant Regge parameter is 2mpv/{Q'^ + m^) 
rather than the 1/x. Consequently, in the small-Q^ region the Regge parameter Xq/x in 
eqs.(0) and (§) must be substituted by {xq/x){1 + rn?p/Q'^) . 

One more remark on kinematics is in order. The BFKL-Regge expansion (^) holds at 
small X < 10~^. In order to model the sea contribution at larger x we multiply (H) by the 
familiar factor (1 — x)"*, with m = 5. This factor does not affect the diffraction region but 
strongly suppresses production of gluons with x > 0.1. 

In applications it is convenient to work with /n(Q^) represented in an analytical form. 
For the leading singularity we have 



1 + RlQ^ 



1 + Co log(l + rlQ^ 



70 



(9) 



which has the large-Q asymptotics |T6 



/o(Q^) (X [as{Q')\ 



-70 



7o 



3Ao 



(10) 



For n > 1 the functions fniQ"^) can be approximated by 

1 _|_ V2f^2 rimax ( 



where 



1 + c„ log(l + r„Q ) - 1, 7„ = 7o5„ 



7n 



(12) 



and n^ax =min{n, 2}. 

Since the relevant variable is a power of the inverse gauge coupling the nodes of fniQ"^) 
are spaced by 2-3 orders of magnitude in Q^-scale and only the first two of them are in 



3 



the accessible range of [|n|, [T^. The first nodes of sub-leading fniQ"^) are located at 



g2 ^ 20 - mCeV^, the second nodes of /slQ^) and h{Q'^) are at ^ 5 • W GeV^ and 
— 2 ■ lO^GeV"^, respectively. The parameterization tuned to reproduce the numerical 
results for fniQ"^) at < 10^ GeV^ is given by eg. ([TT|) . For n = 3 we take a simplified form 
with only two first nodes, because the third node of fziQ"^) is at ~ 2 ■ 10^ GeV'^, way beyond 
the reach of accelerator experiments at small x. The found parameters are listed in the Table. 



n 


On 




rl , GeV-2 


Rl , GeV-2 


.(1) 




5n 


A„ 





0.0232 


0.3261 


1.1204 


2.6018 








0.40 


1 


0.279 


0.1113 


0.8755 


3.4648 


2.4773 




1.0915 


0.220 


2 


0.195 


0.0833 


1.5682 


3.4824 


1.7706 


12.991 


1.2450 


0.148 


3 


0.471 


0.0653 


3.9567 


2.7756 


1.4963 


6.9160 


1.2284 


0.111 



Asymptotically, at 1/a; — > oo, the expansion (||) is dominated by the term fo{Q^){xo/x)^° . 
At moderately small x the sub-leading terms are equally important since A„ ~ l/n. How- 
ever, as it has been pointed out in |1TT| , |l^], for < 10^ GeV^ all /n(Q^) with n > 3 are 



very close in shape to each other. Then we arrive at the truncated expansion 



F2{x,Q' 



3 

E 

n=0 



(13) 



where the term f3{Q'^){xo/x)^^ with the properly adjusted weight factor, a^, stands for all 
terms with n > 3. The addition of this "background" term in eq.(|13]) improves significantly 
the agreement with data for large thus expanding the applicability region of eq. (|13D over 
the whole small-x kinematical domain of HERA. 

The need for a soft pomeron contribution ^2^°-^* in addition to the perturbative gBFKL 
SF's described previously is brought about by phenomenological considerations. A viable 
gBFKL phenomenology of the rising component of the proton structure function over the 
whole range of studied at HERA (real photo-absorption included) is obtained if one starts 
with the Born dipole cross section ^"^(r) as a boundary condition for the gBFKL evolution 
at Xq = 0.03 [|l3, |l^. However, such a purely perturbative input, (^^(r), with Rc = 0.27 fm 
strongly underestimates the cross sections of soft processes and the proton SF at moderate 
Q"^ ~ IGeV^. Therefore, at r > -Rc, the above described perturbative gBFKL dipole cross 
section apt{x,r), must be complemented by the contribution from the non-perturbative soft 
pomeron, anpt{x, r). In terms of the relationship [|l^ between a{x, r) and the gluon structure 
function of the proton, G{x, Q"^), the non-perturbative dipole cross section crnpt{r) at r > i?c 
must be associated with soft non-perturbative gluons in the conventional G{x,Q'^). The 
contribution to G{x, Q^) from the non-perturbative transverse momenta k'^ < Qo ^ rn?p 
persists at all and must not be subjected to the DGLAP evolution. 

Because the BFKL rise of (j(x, r) is due to production of s-channel perturbative gluons, 
which does not contribute to o-nptij) in |1^, |T^ we argued that to a first approximation 
one must consider the energy independent <Jnpt{r) and additivity of scattering amplitudes 
from both the hard BFKL and soft non-perturbative mechanisms. For recent suggestions 
to identify our anptir) with the soft pomeron of the two-pomeron picture see |]18|, [T^. In 



the models of soft scattering via polarization of the non-perturbative QCD vacuum PO, pi 
o'nptir) is interpreted in terms of the non-perturbative gluon distributions. 
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To our opinion, the recently encountered troubles with the small-Q^ extrapolations of 
DGLAP evolution |[ and the failure of DGLAP fits in the Caldwell plot 0, [g are due to 
illegitimate enforcing the DGLAP evolution upon the non-perturbative glue. 

The non-perturbative term F2°^^{Q'^) in eq.(p!3D calculated from eq.(^ ) with a = anptir) 



from [22 1 can be parameterized as follows 



P^"''^Q')= \Ir2q2 [l+clog(l + r^Q^)] , (14) 

where h = 0.1077, c = 0.0673, = 6.6447 GeV'^ and = 7.0332 GeV-^. So, its logQ^. 
derivative levels off at very small ~ 0.15 GeV^ and does not contribute to the observed 
growth of dF2 / d log . 

In Fig.l we confront our estimates to both the HERA data and the fixed target data. In 
the logarithmic slope, dF2/d\ogQ'^ , is derived from data by fitting F2 = a + blogQ^ in 
bins of fixed x. The average value of Q^, (Q^), in each x-bin is derived from the F2 weighted 
mean value of log in that bin. 

As we have noticed above, at moderately small x ~ 10~^ — 10~^ the contribution of the 
sub-leading poles to -^2(0;, Q^) is still substantial (the relative weight factors, a^, with n > 1 
are presented in the table), but toward the region of x ~ 10~^ the leading pole contribution 
starts to prevail. At small the ratio of log Q^-derivatives, r„ = fn/fo, can be estimated as 



Al / X \^o-A„ 



(15) 



where 



Al ^ l/{Rl - Wo') . (16) 

Because the sub-leading SF's, fniQ"^), have node at ~ 20 — 60 GeV^ |Tl|, |12l, their 
contribution to the slope dF2/d\ogQ^ vanishes at Q' ~ 5 — lOGeV^, which is very close 
to the turn-over point in the HERA data. Hence dF2/d\ogQ^ at small follows closely 
9/0/9 log Q2. From (D it follows that at small Q^^ fo{Q'^) behaves like ~ QV(Ag + g2) ^ith 
Aq ^ 0.72 GeV^ coming from (|I6|). Therefore, dF2/d\ogQ'^ rises with up to Q"^ ~ IGeV^ 
then levels off. Only at large Q', when the sub-leading terms enter the game, 9^2/9 log Q' 
decreases and even becomes negative valued at large x. Our estimates shown in Fig. la 
are in good agreement with HERA data 0. The curves are somewhat wiggly because the 
X — (Q^) correlation of the experimental data is non-monotonous one. 

In Fig. lb we compare our predictions with the fixed target data |]l|. Variation of the 
slope in this case is less pronounced since the starting value of (Q^) is {Q"^) ~ 0.54 GeV^ at 
X ~ 10"^ (compare with Q"^ = 0.12 GeV^ at x = 2.1 ■ 10"^ at HERA). It can easily be seen 
that the derivative dfo/dlogQ^ at such is a rather slow function of Q^. The agreement 
of our estimates with the fixed target data is quite reasonable, though there is a systematic 
discrepancy at small x. We recall that there is a certain mismatch between the E665 and 
Hl/ZEUS data on F2{x,Q'^) in the close (x,^^) bins (see Fig.2a and Fig.2b). 
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Figure Captions 



Fig.l Caldwell's plot of 9^2/9 log for the ZEUS data (Fig. la) and fixed target data 
[H (Fig. lb). Our predictions (BFKL-Regge) are shown by the solid lines. Shown by 
the dashed lines is the leading BFKL pole approximation (LPA). 

Fig.2 Description of the HI, ZEUS and E665 F2(x, Q^) data by the BFKL-Regge expansion 
(1): the large-Q2 data (Q^ =3.5, 12, 25, 65, 120 and 200 GeV^) are shown in Fig.2a, 
the small-Q2 data (Q^ =0.11, 0.20, 0.40, 0.65, 0.85 and 1.2 GeV^) are in Fig.2b. For 
display purposes we have multiplied F2 by the numbers shown in brackets. 
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